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Abstract 

The nucleon-nucleon interaction is investigated by using the ImQMD model with the three 
sets of parameters IQ1, IQ2 and IQ3 in which the corresponding incompressibility coefficients of 
nuclear matter are different. The charge distribution of fragments for various reaction systems are 
calculated at different incident energies. The parameters strongly affect the charge distribution 
below the threshold energy of nuclear multifragmentation. The fragment multiplicity spectrum for 
238 U+ i97 Au at 15 AMeV and the charge distribution for 129 Xe+ 120 Sn at 32 and 45 AMeV, and 
197 Au+ 197 Au at 35 AMeV are reproduced by the ImQMD model with the set of parameter IQ3. It 
is concluded that charge distribution of the fragments and the fragment multiplicity spectrum are 
good observables for studying N-N interaction, the Fermi energy region is a sensitive energy region 
to explore the N-N interaction, and IQ3 is a suitable set of parameters for the ImQMD model. 
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I. INTRODUCTION 



The nucleon-nucleon (N-N) interaction is the most fundamental problem in nuclear 
physics. It is related to many nuclear properties and nuclear reactions mechanism, for ex- 
ample, the binding energies, incompressibility coefficient, nuclear structure, fusion-fission 
reactions and so on. Therefore, knowing the N-N interaction potential form is of great 
significance for us to explain the nature of nucleus and study nuclear reaction mechanism. 
The Skyrme force is an effective N-N interaction with various parameterizations that have 
been proposed to the G-matrix for nuclear Hartree-Fock calculations reproducing the basic 
nuclear structure. Since the first work of Vautherin and Brinkfl], who performed fully mi- 
croscopic self-consistent mean-field Hartree-Fock calculations with the Skyrme type effective 
nucleon-nucleon interaction [2-4], many different parameterizations of the Skyrme interac- 
tion have been proposed and have been reproduced the data of the nuclear masses, radii and 
other physical quantities, such as SkP, SkM*, Slyl-7[5-7]. The improved quantum molecular 
dynamic (ImQMD) model adopts the Skyrme type effective interaction and is successfully 
used for intermediate-energy heavy-ion collisions and heavy-ion reactions at energies near 
the Coulomb barrier [8-11]. By combining the known Skyrme forces and the experimental 
data for fusion reactions and heavy-ion collisions at intermediate energies, three different 
parameter sets IQ1, IQ2 and IQ3 were proposed[9, 12, 13]. The different ImQMD parame- 
ters mean different N-N interaction, and they are associated with different nuclear equations 
of state (EOS). 

In low energy heavy-ion collision, the reaction system has a low excitation energy, one 
observes the emission of light particles plus an evaporation residue for light systems or fis- 
sion caused by Coulomb repulsion for heavy systems. As the incident energy increases to 
a threshold energy, the excitation energy of the reaction system will reach maximum limit. 
It implies that the nuclear reaction will enter the multifragmentation process. In order to 
investigate N-N interaction, we will study the charge distribution of fragments for various 
reaction system at different incident energies and the fragment multiplicity spectrum for 
238 U+ 197 Au at 15 AMeV with the ImQMD model by adopting different parameters. The 
structure of this paper is as follows: In Sec. II, relation of EOS and wave-packet width is 
briefly introduced. In Sec. Ill, three sets of ImQMD parameters IQ1, IQ2 and IQ3 will be 
employed to calculate charge distribution for 40 Ca+ 40 Ca with incident energy from 10 to 45 
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AMeV, 129 Xe+ 120 Sn at 32 and 45 AMeV, 197 Au+ 197 Au at 35 and 50 AMeV, and fragment 
multiplicity spectrumat at 15 AMeV for 238 U+ 197 Au. Finally, a summary is given in Sec. 
IV. 

II. EOS AND WAVE-PACKET WIDTH 

The equation of state (EOS) of nuclear matter plays an important role in studying of 
nuclear properties, heavy-ion collisions, neutron stars and supernova. The EOS depends on 
the interactions and properties of the particles in the matter. It describes how the state of 
the matter changes under different conditions. For cold nuclear matter, the EOS is usually 
defined as the binding energy per nucleon as a function E/A = e(p,T,5) of the density p, 
temperature T and and isospin asymmetry 5 = (p n — p p ) / p. For the symmetric nuclear 
matter 5 = [14]. Empirical values of the nuclear matter EOS such as energy per nucleon, 
incompressibility and saturation density of symmetric nuclear matter at T = MeV suggest 
that Eo ~ —16 MeV, A^ « 230 MeV and po ~ 0.16 fm -3 . However, the uncertainty of 
nuclear equation of state still causes some difficulties for an unambiguous determination of 
the model parameters. In the ImQMD model, the equation of state of the symmetric nuclear 
matter can be expressed as [12, 15] 

E (p) c 2/3 . 1 P . a P 1 , P v m 
— j- = ic k p 1 + -a— + p- — -—^ + g T — . 1 

A 2 p (7+IJP0 Po 

Where £ = c /ck, = | ^(^f-) 2 ^ 3 = 75.0MeV fm~ 2 . The coefficient Co can be determined 
by the kinetic energies of nuclei at their ground state [15]. The density distribution function 
p of a system can be read 

Fig. 1 shows the energy per nucleon of symmetric nuclear matter as a function of p/po- It is 
seen that the EOS with the IQ3 is clearly harder than that with the IQ1 and IQ2 in region 
p/po > 1 with the increase of the density. For the parameter set IQ3, the incompressibility 
coefficient is 226 MeV and much higher than the corresponding values for IQ1 (165 MeV) 
and IQ2 (195 MeV). In heavy-ion collision process, the projectile and target firstly contacted, 
compressed and then expanded. The soft nuclear matter have a higher compression density 
and more violent expansion process. While for the finite nuclear system, the expansion 



3 



I I I I I I I I I I I 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 




FIG. 1: The energy per nucleon of symmetric nuclear matter for different parameters. 

process will be influenced by the nuclear surface energy, Coulomb energy and wave-packet 
width, etc. 

In the ImQMD model, each nucleon is described by a coherent state of Gaussian wave 
packet. The system-size-dependent wave-packet width in coordinate space is given by the 
formula: 

a n r = a + ^A 1 / 3 , n = {p,t}. (3) 

Here, a? (a*) denotes the wave-packet width for the nucleons which belong to the projectile 
(target). A p and A t denote the mass number of the projectile and target, respectively. For 
symmetric reaction systems, a r = a% = <r*. The parameter <jq and o\ have been listed in 
Table 1 . The wave-packet width is useful for exploring the influence of the interaction range 
of nucleons and the finite-size effect of nuclei. In Fig. 2, we show the wave-packet width of 
nucleon as a function of system-size with the three sets of parameters. One can see from 
the figure that there is a large difference with the different system. The wave-packet width 
given by the IQ1 and IQ2 obviously is higher than the IQ3 with increasing mass number. 
The three sets of parameters are listed in Table 1. 
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FIG. 2: The wave-packet width as a function of system-size with the three sets of parameters. 



TABLE I: The model parameters. 





a (3 


7 


go 


9t 


V 




Po 


Co f 1 




(MeV) {MeV) 




(MeVfm 2 ) 


(MeV) 




(MeV) (fm 2 ) 


(fm' 3 ) 


(fm) (fm) 


IQ1 


-310 258 


7/6 


19.8 


9.5 


2/3 


32.0 0.08 


0.165 


0.49 0.16 


IQ2 


-356 303 


7/6 


7.0 


12.5 


2/3 


32.0 0.08 


0.165 


0.88 0.09 


IQ3 


-207 138 


7/6 


18.0 


14.0 


5/3 


32.0 0.08 


0.165 


0.94 0.018 



III. RESULTS 

Based on the ImQMD model, the charge distribution of fragments for different reaction 
systems will be employed to explore the N-N interaction. We present the charge distribution 
of fragments for 40 Ca+ 40 Ca at incident energy of 35 AMeV in Fig. 3. The solid triangles, 
solid squares and solid stars denote the results with the IQ1, IQ2 and IQ3, respectively. 
The open circles denote the experimental data[17]. Here we create 500 events for central 
collisions and for each event we simulate the whole collision process until t=3000 fm/c 
without combining statistical model. From the Fig. 3a-3c, one sees that the experimental 
data can be reproduced very well with the ImQMD calculations with three sets of parameters. 
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FIG. 3: The charge distribution of fragments from three sets of parameters for Ca+ Ca at 
incident energy of 35 AMeV. 

However, their nuclear equations of state and the wave-packet width are different. Why do 
they have a similar charge distribution? 

For further investigating the effect of N-N interaction on the heavy-ion collisions, we 
studied the charge distribution of fragments for different reaction system at different incident 
energies. In Fig. 4 we show the charge distribution of fragments calculated by the ImQMD 
model with three sets of parameters for 40 Ca+ 40 Ca at incident energies from 10 to 45 AMeV. 
One sees that the results from IQ1, IQ2 and IQ3 for the charge distribution of fragments are 
relatively close to each other at the energy region 35-45 AMeV, but there exist significant 
differences at the energy region 10-30 AMeV. One can also see that all of the peak from the 
IQ1, IQ2 and IQ3 become gradually lower with increasing incident energies at the energy 
region 10-30 AMeV. In this region, it is the transition region of the fusion reaction to 
the multifragmentation for 40 Ca+ 40 Ca reaction system. The fusion evaporation of reaction 
system is strongly affected by the ImQMD parameters. Above 30 AMeV, the nuclear reaction 
turn to the multifragmentation, the distribution of fragments is not sensitive to the ImQMD 
parameters. 

Fig. 5 shows the charge distribution of fragments for 129 Xe+ 120 Sn at 32 and 45 AMeV 
[18], 197 Au+ 197 Au at 35 and 50 AMeV[19], which are calculated by the ImQMD model 
using three sets of parameters. One can see that the charge distribution of fragments at 45 
AMeV for 129 Xe+ 120 Sn and 50 AMeV for 197 Au+ 197 Au are also relatively close to each other 
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FIG. 4: The charge distribution of fragments from three sets of parameters for 40 Ca+ Ca at 
incident energies from 10 to 45 AMeV. The solid circles, open circles and open triangles denote 
the results with the IQ1, IQ2 and IQ3, respectively. 

calculated by three sets of parameters. While there still exist significant differences at the 
lower incident energy of 32 AMeV for 129 Xe+ 120 Sn reaction and 35 AMeV for 197 Au+ 197 Au 
reaction. Compare 129 Xe+ 120 Sn and 197 Au+ 197 Au reaction system with Ca+ Ca, we 
find that the sensitive energy region for the ImQMD parameters is depend on the reaction 
systems. The beginning of insensitive energy region for heavy system is higher than the light 
system. It is well known, the mean field plays a dominant role at the low energy region. 
However, the nucleon-nucleon collision becomes more and more important with increasing 
energy. In the transition region, the mean field and the nucleon-nucleon collision work 
together, so it is a sensitive energy region for the N-N interaction. By comparing with the 



7 



10' 
10 1 



o 129 Xe+ 120 Sn, 32 AMeV 




10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 80 



b i97 Au+ i97 Au 35 AMeV 




o exp 

a IQ1 

■ IQ2 

* IQ3 



N 
T3 



10 1 
10° 
10 1 
10" 2 
10 3 
10 4 
10 s 



129 Xe+ 120 Sn, 45 AMeV 




an o 



i 1 1 1 1 ■ 1 - 



10 20 30 40 50 




10 20 30 40 50 



FIG. 5: The charge distribution of fragments for 129 Xe+ 120 Sn at 32 and 45 AMeV, 197 Au+ 197 Au at 
35 and 50 AMeV. The solid triangles, solid squares, open stars and open circles denote the results 
with the IQ1, IQ2, IQ3 and the experimental data, respectively. The 197 Au+ 197 Au at 50 AMeV 
does not give the experimental data. 

experimental results, we find that the the experimental data can be always reproduced very 
well by the IQ3. The IQ3 is a suitable set of parameters in the ImQMD model for heavy- ion 
collisions. 

In order to further test the above conclusion, we calculated fragment multiplicity spec- 
trum of 238 U+ 197 Au reaction by using three sets of parameters. In Fig. 6, we show the 
comparison of calculation results with the experimental data for fragment multiplicity spec- 
trum of 238 U+ 197 Au reaction at 15 AMeV. The open squares, open triangles, solid stars and 
open circles denote the results with IQ1, IQ2, IQ3 and the experimental data[20], respec- 
tively. We created 4000 events from central to peripheral collisions and counted the number 
of fragments in each reaction except in the case where the charge Z < 8 as that did in the 
experiment. The experimental data shows that two-body events exhaust only about 5% and 
the three- and four-body events exhaust approximately 83% of the total reaction events. 
From Fig. 6, one can see that the ImQMD calculation result reproduce the experimental 
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FIG. 6: The fragment multiplicity spectrum of U+ Au reaction at 15 AMeV. 

multiplicity spectrum reasonably well with the IQ3 rather than the IQ1 and the IQ2. It 
implies that the IQ3 is reasonable in the ImQMD model for heavy-ion collisions. 



The N-N interaction has been investigated by using the ImQMD model with the three 
sets of parameters IQ1, IQ2 and IQ3. The charge distribution of fragments at different 
incident energies were calculated based on the three sets of parameters. The calculation 
results demonstrate that the charge distribution and the fragment multiplicity spectrum are 
good observables and the Fermi energy region is a sensitive energy region for probing the 
N-N interaction. The charge distribution of fragments below the threshold energy of nuclear 
multifragmentation are very sensitive to the ImQMD parameters. By comparing the charge 
distribution of fragments and fragment multiplicity spectrum with the experimental data, 
we find that the IQ3 is a suitable set of parameters for the ImQMD model. 
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